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Abstract 

Multiplicities of periodic orbit lengths for non-arithmetic Hecke 
triangle groups are discussed. It is demonstrated both numerically 
and analytically that at least for certain groups the mean multiplicity 
of periodic orbits with exactly the same length increases exponentially 
with the length. The main ingredient used is the construction of joint 
distribution of periodic orbits when group matrices are transformed 
by field isomorphisms. The method can be generalized to other groups 
for which traces of group matrices are integers of an algebraic field of 
finite degree. 

1 Introduction 

For chaotic systems the density of classical periodic orbits with a given length 
increases exponentially. In particular, for all constant negative curvature 
surfaces generated by discrete groups one has the universal asymptotics (see 

e.g. ini) 

Ptotai(0 ^ y • (1) 

Much less is known about multiplicities of periodic orbits with exactly the 
same length. Usually it is assumed that the mean length multiplicity of 
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periodic orbits for generic systems depends only on exact symmetries and for 
models without geometrical symmetries the mean multiplicity g equals 2 or 
1 for systems respectively with or without time-reversal invariance. 

Physically it means that, in general, there exists no reason that two dif- 
ferent periodic orbits would have the same length except for time-reversal 
invariant systems where almost all trajectories can be traversed in two oppo- 
site directions which implies that g = 2. In semiclassical approach to spectral 
statistics of chaotic systems the distinction between these two classes of mod- 
els is reflected in different behaviour of the two-point correlation form factor 
at the origin which agrees with the predictions of the random matrix theory 

HI, 0. 

For the free motion on constant negative curvature surfaces generated by 
discrete groups the situation is different. In such hyperbolic models periodic 
orbits are in one-to-one correspondence with conjugacy classes of group ma- 
trices and the length of a periodic orbit, Ip, is directly related with the trace 
of a matrix M representing each class (see e.g. ^T]) 

iTrA/l-/ 2cosh/,/2, ifdetM = l 

I ™ '"12 sinh lp/2 , if det M = -1 " 

Hence, any relations between traces of group matrices imply connections 
between periodic orbit lengths. 

The extreme case corresponds to the so-called arithmetic groups (see e.g. 
and references therein). For such groups traces of group matrices can take 
only quite restricted set of values and the number of possible traces less than 
a given value is asymptotically |lj 

N{\TrM\ < X)^-^ CX (3) 

with a system dependent constant C. Because X '-^^ e'^^ but not all possible 
values of traces really appear for group matrices, the number of periodic 
orbits with different lengths when Ip ^ oo has the following upper bound jl] 

iVdiff.(/p<0 <Ce'/^ (4) 

Define the mean multiplicity of periodic orbit length as the ratio of the density 
of all periodic orbit to the density of periodic orbits with different lengths 

Ptotal(0 fr\ 

^^^^ = ^wn 
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where Pdifr.(0 = dN^iff {lp < l)/dl. 

From the above formulas one proves |3] that for arithmetic groups the 
mean multiphcity is exponentially large and has the following estimate from 
below 

m) > ^ . (6) 

In classical mechanics such large multiplicities play a minor role but their 
interference changes drastically quantum mechanics of arithmetic groups. In 
particular, spectral statistics of arithmetic systems is close to the Poisson 
statistics typical for integrable systems and not to the random matrix statis- 
tics conjectured for chaotic models jl], 0, [SI- 

Arithmetic systems are very exceptional but according to the Horowitz- 
Randol theorem ^2); Q2| for all hyperbolic models generated by discrete 
groups multiplicities are unbounded. Nevertheless, multiplicities covered by 
this theorem are quite rare and a priori assumption would be that for non- 
arithmetic hyperbolic models the mean multiplicity equals 2 as for generic 
time-reversal invariant systems. 

Numerical calculations performed in P| indicated that it is not always the 
case. In that paper certain non-arithmetic Hecke triangles were considered 
and it was observed that mean multiplicity of periodic orbits with length I 
seems to increase exponentially 

-g{l) ~ e^' (7) 

with an exponent A < 1/2. 

The purpose of this paper is two-fold. First, we perform numerical calcu- 
lations of periodic orbits for much larger lengths that in and, second, we 
develop a method which gives a lower bound of multiplicities, thus in certain 
cases confirming analytically exponential growth of multiplicities. 

The plan of the paper is the following. In Section |21 we discuss general 
properties of Hecke triangle group matrices. In Section |21 results of numerical 
calculations of periodic orbit length multiplicities for a few Hecke triangles 
are presented. As traces of Hecke triangle group matrices are integers of an 
algebraic field, each group matrix defines not one but a few different lengths 
corresponding to different isomorphisms of the basis field. In Section E] the 
construction of the joint distribution for periodic orbits with all transformed 
lengths fixed is discussed. In Section |S]it is demonstrated how the knowledge 
of this joint distribution permits to calculate the lower bound of the periodic 
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orbit length multiplicity. In Section IF!T] the computations are performed for 
the simplest case of Hecke groups with n = 5 , 8, 10, 12 which are charac- 
terized by the existence of only one non-trivial isomorphism. In Appendix A 
it is proved that for Hecke triangle groups all transformed lengths are smaller 
than the true length. This inequality is sufficient to ensure that for Hecke 
groups with only one non-trivial isomorphism length multiplicities increases 
exponentially. Our results agree well with direct numerical computations of 
periodic orbit multiplicity for these groups. In Section other Hecke groups 
are shortly considered. It appears that in all investigated cases except groups 
with one non-trivial isomorphism length multiplicities increase so slowly that 
the direct check is practically impossible. In Section El we briefly discuss the 
influence of periodic orbit length multiplicities on the spectral statistics for 
corresponding systems. In Section [7| a summary of the results is given. In 
Appendix B a saddle point method of calculation of the joint distribution of 
periodic orbit lengths is discussed. 



2 Hecke triangles 

Hecke triangles are hyperbolic triangles with angles 0, 7r/2, vr/n with integer 
n > 3. All of them are fundamental regions of discrete groups G„ generated 
by reflections across its sides. Let us denote the reflection across the side 
connecting angles and 7i/2 by A, the one across the side connecting angles 
and tt/u by B and the last one by C. From geometrical considerations 
these transformations obey the deflning relations 

= B^ = C^ = 1, {ACf = {BCY = -1 . (8) 

The explicit form of A, B, and C can be chosen as follows 

).c^(;j). (0, 

Arbitrary matrix from the Hecke triangle group Gn is a word of these letters. 
Due to (jHl) these symbols have a complicated grammar. For our purposes it 
is convenient to introduce new symbols 

a™ = C{ABr , 
= CB{ABr , 
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a™ = CBCiAB)"" , (10) 



a 



n-2 



CBCB...{ABy 



n—2 symbols 



where m = 1, 2, . . . are positive integers. 
Explicitly up to unessential overall sign 

m ( -«fc amak + ak-i\ „ f -ak-i amak-i + at \ 

^ V ~^k+i amak+i + at ) \ -a^ amat + flfc+i / 

where from now we denote 

a = 2 cosvr/n (12) 
and Ofc = afc(Q;) are the Chebyshev polynomials of the second kind of a 

sin(fc7r/n) 
sm(7r/n) 

Using the defining relations (jHI) one can proves j3] that conjugacy classes in 
Gn (and, consequently, periodic orbits in Hecke triangles) can be constructed 
as free words in these new symbols with the only restriction that cyclic per- 
mutations correspond to the same orbit. 

Due to specific form of generators matrix elements of the Hecke trian- 
gular group matrices are polynomials with integer coefficients of the variable 
a = 2cos7r/?T,, thus forming naturally a subfield of the cyclotomic field of 
degree 2n. 

The constant a defined in (fT^ obeys a polynomial equation Pat (a) = 
with integer coefficients 

Pn{x)= n (x-afc(n))=x^ + ... (14) 
k = odd 

{k,n) = 1 

where 

ttfc(^) = 2cos{TTk/n) (15) 

and the product is taken over all odd integer k coprime with n. The total 
number of such integers and, consequently, the degree of the defining equation 
is 

N = ^ip{2n) (16) 
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where <f{p) is the Euler (^-function which counts the number of integers less 
than p and coprime with p. 

In Table ^ the explicit forms of the defining polynomials for low values 
of n are presented. In the last column of this table we give for later use the 
discriminant of these polynomials defined as the square of the product of all 
roots 

A„ = Yl [ak{n) - am{n)f (17) 

k<m<N 

where ak{n) is given by (fT3j) . and the product is taken over all odd integers 
k < m both coprime with n. For even n and odd k an-k{.n) = —ak{n), and 

A„ = 2^A(^)Ai°) (18) 

where A^^'"^ are the discriminants of even and odd powers of a 

k<m<N/2 

Ai°^ = n [Mn)am{n){al{n)-al{n))f . (20) 

k<m<N/2 

Therefore all matrix elements and, in particular, traces of Hecke group 
matrices have the following form 

N-l 

TrM = J2 ^ka' (21) 

A:=0 

with integer coefficients Uk- 

As matrix elements of Hecke groups are algebraic integers of the totally 
real field (Q(2cos7r/n) it is natural to consider in parallel all isomorphisms 
of this field defined by the following substitutions 

TT A/ 

(fk : a — > ak = 2 cos — (22) 

n 

for all odd integers k < n coprime with n. 

In general, the number of such isomorphisms equals the degree of the 
defining polynomial but in our case a and —a both correspond to the same 
group. Hence, when this transformation belongs to the group of isomor- 
phisms (which is the case for even n), it does not change periodic orbit 
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Table 1: Irreducible monic polynomials defining the field Q(2cos7r/n) for small 
n. N is the degree of the polynomial, A„ is its discriminant. For n even the 
discriminant is given as the product of 3 terms. The second and third factors 
represent discriminants for even and odd functions (see H19|) and ()2U() ). 



n 




Pn{x) 


A 


5 


2 


x"^ — X — 1 




5 




7 


3 


— — 2x + 1 




72 




8 


4 


x^ -4x2 + 2 


24 


■23. 


24 


9 


3 


x^ — 3x — 1 




34 




10 


4 


x^ — Sx^ + 5 


24 


• 5 ■ 


52 


11 


5 


x^ - x'' - 4x3 + 3x2 + 3x - 1 




114 




12 


4 


x^ -4x2 + 1 


24 


■ 12- 


12 


13 


6 


x^ — x^ — Sx^ + 4x3 + 6x2 _ _ I 




134 




14 


6 


x^ - 7x4 ^ ^4^2 _ Y 


26 


■72. 


73 


15 


4 


x^ + x3 — 4x2 — 4x + 1 




3253 





lengths. Consequently, the dimension of the group of isomorphisms of peri- 
odic orbit lengths, q, is 

^ i N, for odd n , . 

^ 1 |A^, for even n ' ' 

In particular, the following four cases corresponds to the simplest case of the 
groups of isomorphisms of the order 2 (cf. Table 

n = 5 , 8 , 10 , 12 . (24) 

It appears that multiplicities of periodic orbit lengths depend strongly on the 
number of isomorphisms so we consider first the case 



3 Numerical Calculations 

At Fig. ^ we present the numerically computed multiplicity for the Hecke 
triangles (0, 7r/2, vr/ri) with n = 5,8,10,12 for lengths / < 20. White lines 
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Figure 1: Mean multiplicities of periodic orbit lengths for Hecke triangles 
(0, 7r/2, vr/n) with (from top to bottom) n = 12, n = 5, n = 8, and n = 10 
for I < 20. White hues are numerical fits 



indicate a two-parameter fit to these data in the form g{l) ^ a„e 



bnl 



n = 


5 : 


m) 


^ 1.235e°-ii^' 


n = 


8 : 


m) 


^ 1.095e°-ii^' 


n = 


10 : 


m) 


^ 1.143e°-°6^' 


n = 


12 : 


m 


^ 0.986e°-^^°' 



(25) 
(26) 
(27) 
(28) 



Expressions ()25 p - ()28|) fit numerical data in the given interval of lengths 
pretty well. But they are purely best least-square numerical fits and no 
attempts were made to determine the accuracy of coefficients. In Section EH] 
it is demonstrated that our approach suggests different formulas for these 
quantities (see (f75)) ) which, nevertheless, are practically indistinguishable 
from the above simple expressions in the considered interval of lengths (cf. 
Fig. ED. 

For larger lengths the exponential proliferation of periodic orbits makes it 
difficult to compute and store in the memory all periodic orbits. Nevertheless 
the determination of periodic orbits in a reasonably short interval of lengths 
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Figure 2: Mean multiplicities of periodic orbit lengths for the Hecke triangle with 
angles (0, 7r/2, 7r/5) for / < 25. Solid line is the fit (j^ . 

is still feasible. At Fig. |21we present the result of the numerical computation 
of the length multiplicity for the Hecke triangle with n = 5 till / = 25. Each 
small circle at this figure for / > 20 corresponds to one million of periodic 
orbits. The solid line is the fit ()25|) obtained from data at small /. It is clearly 
seen that accuracy of the fit does not change noticeably with the increasing 
of periodic orbit lengths. 

4 Length distribution for different 
isomorphisms 

For the Hecke triangle groups (and for certain other groups as well) traces of 
group matrices are integers of an algebraic field of finite degree. Therefore 
each group matrix M gives rise not only to one usual length Q but to q 
different lengths corresponded to q different isomorphisms of the basis field 
applied to a matrix M. Asymptotically 

/fc = 21n|Tr(^fc(M)| . (29) 
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In this definition li corresponding to the identity transformation is the true 
length of a periodic orbit and all other with k >2 are additional quantities 
which we call transformed lengths. 

For arithmetic systems (see e.g. transformed traces are restricted 

|Try;fc(M)|<2 (30) 

for all A; > 2 which leads to the very large length multiplicities lengths for 
such groups (jni). 

The main ingredient of our approach to the problem of length multiplicity 
for non-arithmetic groups is the determination of the joint density of periodic 
orbits in intervals lk,lk + dh for all > 1. For clarity we first consider groups 
with only one non-trivial isomorphisms where each hyperbolic group 
matrix permits to define two lengths, h and ^2- 

Let R{li,l2)dhdl2 be the number of periodic orbits with the first length 
in the interval /i, Zi + dh and the second (transformed) length in the interval 
hyh + dl2- Taking into account ((H) one concludes that 

R{hj2)-^P{h,l2) (31) 

where P{li, I2) has the meaning of the probability density of periodic orbits 
with lengths li and I2 normalized such that 

/oo 
P{h,l2)dl2 = l. (32) 
-00 

No general arguments determining the form of P(/i,/2) are known to the 
authors. As li is only one fixed quantity with the dimension of the length, 
from physical considerations it is quite natural to assume that for large li 
and I2 this function has the following scaling form (see also Appendix B for 
another argument) 

P(/i,Z2) =A(Zi)exp[Zi/(/2//i)] (33) 

with a certain (smooth) scaling function f{u) where u is the ratio of two 
lengths. 

When /i ^ 00 the prefactor A{li) can be determined in the saddle point 
approximation from the normalization condition (jH^ 

V27r(T^Zi 
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where cr^ = 1/|/"(Mc)|. Here Uc is the point of the maximum of f{u): f'{u^ = 
0, and f"{uc) is the second derivative of the function f{u) at this point. From 
fl32|) it follows that the value of f{u) at the point of the maximum is zero 

f{uc) = . (35) 

At Fig. El we present numerically computed function P(/i,/2) for the Hecke 
triangle with n = 5 for 10^ orbits near li ^ 19.8 (which corresponds to 
|Tr M\ = 20000) together with the Gaussian fit to the data in the form 

P{h,k) = aoexp • (36) 

The least square fit gives the following values of the parameters 

ao = .116 , A = 6.69 , = 11.74 . (37) 

At Fig. EJthe best fit values of A(/i) and cr^(/i) are given for different values 
of li- The data are linear on li and can be approximated by the following 
straight lines 

A = .330/i + .187 , = M6h - -337 (38) 

which support the scaling ansatz 

The peaks at Fig. IHl correspond to words in the code ()11|1 with a small 
number of letters a™ but with big values of m's. Ignoring all elements except 
the ones multiplied by the largest possible numbers of m's one can approxi- 
mate the periodic orbit length as follows 

/ ^ 2 ln(a^ mi . . .rup . . . a^p) • (39) 

Hence, in this approximation the difference between transformed lengths and 
the true length is a finite constant 

/fc - /i ^ 2 ln(v9fc(aP Ofc, . . . CfcJ) - 2 ln(a*' Uk, . . . a^J . (40) 

To compute the joint distribution of transformed lengths one considers pe- 
riodic orbits with the true length confined in a small interval. The above 
expression means that orbits corresponding to small numbers of initial sym- 
bols have transformed lengths at finite distances from /i and, consequently, 
they produce peaks at these distances. The quality of such approximation 
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Figure 3: Probability density of transformed length with li = 19.8 for the Hecke 
triangle with n = 5. Thick line represents the Gaussian fit (|36|) . (|37|) . 

quickly deteriorates with increasing of p due to the omitting lower powers of 
m's and in real calculations only a few peaks with small p are visible. 

For the Hecke triangle with n = 5 defined in fjl3p are either 1 or a, 
and all differences between two lengths are 



with integer m which agree well with the positions of the peaks at Fig. El 

At Figs. IHJIHl we plot numerically computed scaling functions f{u) for 
the Hecke triangles with n = 5,n = 8,n = 10, and n = 12 for different 
intervals of periodic orbit lengths. The curves for different lengths seem to 
be superimposed thus supporting the scaling ansatz (jH^ . Irregular points at 
h/h ~ 1 correspond to the above mentioned peaks pn)) related with words 
with small number of symbols and are irrelevant at large li. 

The scaling functions fih/h) for the Hecke triangles with n = 5 and 
n = 8 are close to each other and can be reasonably well described by the 
following parabolic fit 




(41) 



.094 + .56x - Mx^ 



(42) 
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Figure 4: Gaussian fit parameters (jSHJ versus the length of periodic orbits for the 
Hecke triangle with n = 5. Lower line: A. Upper line: o"^. Solid lines are the linear 
fits (|H5|) to these data. 

indicated by dashed lines at Figs. El and IHl 

The scaling functions for the Hecke triangles with n = 10 and 12 have 
more complicated form. At Fig. U\ the dashed line indicates the cubic fit to 
the data in the interval [.4, 1] 



At Fig. ISlthe dashed line shows the parabolic fit in the interval [.1,1] 



For Hecke triangle groups with n different from (j24j) there exist more than 
one non-trivial isomorphisms and, consequently, the joint distribution of all 
lengths have the form similar to ()31|) but with larger number of transformed 
lengths 




(43) 




(44) 




(45) 



(46) 
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-0.2 0.2 0.4 0.6 0.8 1 

'A 

Figure 5: Scaling function /(h/h) for the Hecke triangle with n = 5. Circles, tri- 
angles, and squares represent data for 10^ orbits near respectively h 19.8 , Zi ^ 
19.47 , li ~ 19.02. Solid line connects points with li ~ 19.8. Dashed line is the 
parabolic fit H42|) to the data with li ~ 19.8. Thick solid line is the straight line 
y = {x-l)/2. 
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l/l, 



Figure 6: The same as at Fig. |21but for the Hecke triangle with n = 8. 
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Figure 7: The same as at Fig. |21but for the Hecke triangle with n = 10. Dashed 
hne is the cubic fit ()43() to the data with li « 19.8 in the interval [.4, 1]. 




Figure 8: The same as at Fig. |21but for the Hecke triangle with n = 12. Dashed 
line is the parabolic fit (|44j) to the data with li ~ 19.8 in the interval [-1,1]. 
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with a certain function f{x2, ■ ■ ■ ,Xg) = f{x) depended only on ratios Xk 
Ik/h and in the saddle point approximation 



det d'^f/dxidxjl 



\ (27r/i)«-i 

where the derivatives are taken at the point of the maximum of f{x) 



(47) 



5 Number of periodic orbits with different 
lengths 

The importance of the knowledge of the joint distribution of periodic orbit 
lengths for all possible isomorphisms is related with the fact that two periodic 
orbits for Hecke triangle groups have exactly the same length iff all their 
transformed lengths are the same. 

Let us consider the simplest Hecke group with n = 5. In this case the 
traces of group matrices, ti and t2 = V^2(^i) can be written as 

ti = no + riiXi , t2 = % + ^1^2 (48) 

where no, tt-i are integers, Ai = 2 cos(7r/5) is an element of our basis field and 
A2 = 2cos(37r/5) is the transformed value of Ai. 

These equations determine the transformation from variables ti, ^2 to vari- 
ables no, ni and 

dtidt2 = Jdnodrii (49) 

where the Jacobian of this transformation is the square root of the discrimi- 
nant (fTTj) of the defining equation 

J= IA2- All = yA^= Vs. (50) 
As ti = e''/^ the precedent equation gives 

driodrii = C5e^'/^+^^/^dhdl2 (51) 

with C5 = l/(4v^). 

Because no and ni are integers this equation means that in a volume dlidl2 
there are at most [Ce'^/^^'^/^j possible values of no, ni {[x] is the integer part 
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of x). This relation signifies that the density of the maximal number of 
periodic orbits with different lengths obeys asymptotically the inequality 

Pdifr.(^i,^2)<C5e'^/'+'^/'. (52) 

We stress that such arguments can give, in principle, the estimate from above 
because not all values of uq and ni are possible for the Hecke group Gn, 
otherwise one obtains the Hilbert modular groups which are discrete groups 
only in higher dimensional complex planes. 

For other Hecke triangle groups with one non-trivial isomorphism (I24j) 
(i.e. for n = 8, 10, 12) the defining equation is of degree 4 but traces of group 
matrices contain either even or odd powers of a 

h = no + n2a^ , or ti = riia + naa^ (53) 

and the result is similar to 

P<^m.{liM)<Cne''"^'''^ (54) 

but with 

Cn = + (55) 

where A^*^'"-* are discriminants ()19|1 and (j2()|) . 
For general Hecke group with q isomorphisms 

p,,,XhJ2,...,Q<Cj^^'^^^/'^---^^^/' (56) 

where for odd n 

Cn = — ^ , (57) 

and for even n 

Cn=^= + ^==- (58) 

Eq. (jl^ means that in a volume dh . . . dig there is 

PtotXn = ^A{h)exp{hf{yiu...JJli)) (59) 
n 
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periodic orbits with all transformed lengths fixed. On the other hand in the 
same volume the maximum number of periodic orbits with different lengths 
is restricted by the inequality 

Pdifr.(r)<ae'^/'+'^/'+-+'^/^ (60) 

Consequently, the maximum number of periodic orbits with different 
lengths is 

(maximum)/, \ _ [ Ji Ji } PdiS.{l ) if Pdifr.(^ ) < Ptot.(^ ) /^^i \ 

Pdiff. lengths IhJ -/ "'2 ••• «fg S ,7'. /Tn^ /Tn " ^^^^ 

J [ Ptot.{l ) if Pdifr.(' ) > Ptot.{l > 

As both densities increase exponentially with li the dominant contribution 
to this integral is given by vicinities of boundary points where 

Pdiff.(^"') =Ptot.(^"') • (62) 

In the leading order of li these points are determined from the equality of 
the exponential factors of these functions 

h + hfih/h, IJh) = \{h + ... + /,). (63) 



Denoting Ik/h by one gets the equation independent of / 



/(x2, . . . , Xg) = i(x2 + . . . + - 1) . (64) 
5.1 Groups with one non-trivial isomorphism 



In the simplest case of groups (j24|) where only one transformed length exists 
Eq. is reduced to the equation of one variable x = X2 

/(x) = ^(x-l). (65) 

In Appendix A it is proved that for the Hecke groups the transformed lengths 
corresponding to all non-trivial isomorphisms are smaller than the true length 

lk<li. (66) 

Consequently, f{x) is situated at the left from the line x = 1 and as /(mc) = 
f{x) is negative when x < 1. As Mc < 1 Eq. ()65|) for groups with one non- 
trivial isomorphism always has a solution x < 1. In Table |21 we present 
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Table 2: Parameters for Hecke triangles with n = 5 , 8, 10, 12. The second 
column is the curvature in the point of the maximum. The third column gives 
the value of C„ in (|54jl . The fourth and the fifth columns are the ordinate of 
the intersection point and the modulus of the slope of /(x) at this point. The 
nest three columns are parameters in (|74|) . The last column gives the numerically 
computed prefactor in ((75|) . 



n 




Cn 










Gn 




5 


.6 


.11 


.74 


.66 


.13 


.35 


1.32 


2.34 


8 


.6 


.15 


.74 


.66 


.13 


.35 


1.11 


2.22 


10 


.55 


.16 


.82 


.82 


.09 


.43 


1.22 


2.46 


12 


.68 


.14 


.64 


.73 


.18 


.39 


1.25 


1.77 



approximate values of this intersection point, Xn, for different values of n 
found from Figs. [SHHl As claimed in all these cases the solution exists and 
Xyi ^ 1 . 

In the next order one can write 

h = xJi + en . (67) 
Expanding Eq. (|62|) to the first order of e one gets 



Cn^J27ialll = exp{-en{kn + .5) + 0{ei/h)) (68) 

where fc„ = |/'(a;„)| is the modulus of the derivatives at the point of the 
intersection. Therefore 



1 



\niCnJ2nallf] . (69) 



kn + .5 

Together these formulas demonstrate that for the Hecke triangles (j^^ at the 
intersection point 

PdifT.(^i, ^2) = Ptot.(/i, k) ^ Dj-^'+^-^'Hf^-'^ (70) 

where 

1 C*^"^" 

^" = 2^;:tt' = (2vra„)'3"/2 • ^^^^ 
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The integration in in the hmit of large li can be performed by parts and 
finally 

, . , 2D ph{l+x„)/2 

(maximum) /? \ '^^n ^ /^n\ 

Pdiff. lengths I'l J = ]^ _ ^3/3„/2 " y'"^) 

The mean multiplicity of periodic orbit lengths is the ratio of the total density 
of periodic orbits to the density of orbits with different lengths. Hence 

m > Gn-^ (73) 

where 

A„ = , z/„ = 1 - , Gn = -^j^ ■ (74) 

At Table El we present approximate values of these parameters computed 
from Figs. EHHl At Fig. IHl we compare data of length multiplicities for the 
Hecke triangles (j^K|l - (f^ with the formula of the form ()73|1 

^(0 = Kn-^ (75) 

with the computed values of A„ and Un from Table El but with a prefactor 
Kn calculated from the best fit to the data (see the last column of Table EJ. 
The 'theoretical' curves ()75|) are practically indistinguishable from the best 
fits Note that fitted prefactors is always bigger than Gn, just 

confirming that estimates (|73p and (j74j) give only lower bounds. Though in 
principle not all integers are allowed in (j2H) . these results seem to indicate 
that in the mean the ratio of the density of allowed integers to all integers 
for the groups is finite. 



5.2 General case 

For general case of g > 2 isomorphisms the arguments, in principle, remain 
the same. One has to perform the following three steps: 

• to check that required solutions of Eq. do exist, 

• to find on {q — 2)-dimensional manifold of these solutions a point with 
the maximum of the sum X2 + ■ ■ ■ + Xq, 

• to compute the integral in a small vicinity of the point of the 
maximum. 
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Figure 9: Comparison of numerically computed fits ()25() - (|28() for length multi- 
plicities for Hecke triangles (dotted lines) with formulas (|75|) with fitted prefactor 
(solid lines). From top to bottom: n = 12, n = 5, n = 8, n = 10. 

For groups with only one non-trivial isomorphism the inequality ()66|1 was 
sufficient to ensure the existence of a solution of Eq. (jMj) . For other groups it 
is not the case and one has to rely mostly on numerical calculations. For ex- 
ample, the necessary condition of the existence of solutions of Eq. ()64|) is that 
at the point U2, ■ ■ ■ , Ug of the maximum of the scaling function /(x2, . . . , Xg) 
the sum U2 + ■ ■ ■ + Ug is less than 1. 

At Fig. Uniwe present the contour plot of the scaling function f{x2,X3) 
for the Hecke triangle group with n = 7 computed from 10^ points near 
h = 25. The contour lines correspond to the sections of the scaling function 
(normalized so that at the maximum it equals zero) at heights —2-10~'^k{2k — 
1) for k = 1, . . . ,9. Numerically from this figure one gets that for the Hecke 
group with n = 7 the solution of equation f{x2,X3) = {x2 + x^ — l)/2 do 
exist and the point with the maximum X2 + x^, corresponds approximately 
to the fourth contour line. It means that the density of maximal number of 
different periodic orbit lengths increases as 




(76) 
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Figure 10: Contour plot of the scaling function /{h/h, h/h) for the Hecke triangle 
with n = 7. Thick line is the solution of equation f{x2,x^) = {x2 + — l)/2. 

where A7 ~ .006. Correspondingly, the mean multiplicity of periodic orbit 
lengths can be estimated (without a prefactor) as 

-g{h) > e-°°^'^ . (77) 

Though it is an exponential increase, the exponent is so small that at really 
accessible lengths h of the order of 20 it practically remains a constant and 
the prefactor dominates. At Fig. ^2 we plot the numerically computed mean 
multiplicity for the Hecke triangle group with n = 7 (averaged over interval 
of traces equal 10). Instead of increasing it shows a slow decrease but the 
best fit to the data in the form 

9 = «y (78) 

gives 

a ^ 3.55 , b ^ .007 , c ^ .168 (79) 

which is larger than ()77|) . Unfortunately, the limited interval of lengths and 
very slow increase of the multiplicity do not permit to obtain clear conclu- 
sions. 
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Figure 11: Mean multiplicity for the Hecke triangle with n = 7. Solid line is the 
fit CZHD, ([7^. 

At Fig. ^1 numerically computed length multiplicities for the Hecke tri- 
angles with n = 9 and n = 11 are presented. Similarly to the n = 7 case 
the data indicate a slow decrease which is more pronounced for the n = 11 
triangle. Is this decrease just a lower-length phenomenon or do multiplicities 
in these cases tend to a constant cannot be answered from the accessible 
data. 

We stress that though the data for the Hecke triangles with n = 7, n = 9, 
and n = 11 do not show clear increase of mean multiplicities they fluctuate 
around values bigger than 2 which differs from the usual expectation. Also 
in all figures we present multiplicities averaged over some length. The true 
multiplicities fluctuate wildly around the mean confirming unusual character 
of the Hecke triangle groups. 

6 Spectral statistics of non-arithmetic Hecke 
triangles 

It is well accepted that length degeneracy of periodic orbits has a profound 
effect on spectral statistics. According to semiclassical theory of spectral 
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Figure 12: Mean multiplicities for the Hecke triangles with n = 9 (the upper 
curve) and n = 11 (the lower curve). White lines represent additional smoothing 
of the curves. 

statistics P], j2] the two-point correlation form factor for chaotic billiards in 
the diagonal approximation is 

K^^'^^^t) = gm)t (80) 

where g{l) is the mean multiplicity of periodic orbits with the length I and 
l{t) = AttH. 

For systems without (resp. with) time-reversal invariance g = 1 (resp. 
g = 2) and (j8(J|) gives the first term of the expansion of the two-point corre- 
lation form factors for standard random matrix ensembles (see e.g. Pj). 

For models considered in the preceding Sections the mean multiplicity 
g{l) increases exponentially as in (f75|) and the form factor calculated in the 
diagonal approximation differs from the random matrix predictions. 

To consider the spectral statistics we compute numerically eigenvalues of 
the Laplace-Beltrami operator with the Dirichlet conditions on the bound- 
aries of the Hecke triangles for different values of n. At Fig. we present 
the differences between the integrated nearest-neighbor distributions and the 
Wigner ansatz for this quantity [Nw[s) = 1 — e"'^ for the Hecke triangles 
with = 5, 7, 12. For comparison on these graphs thick solid lines indicate 
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Figure 13: Differences between integrated nearest-neighbor distributions for the 
Hecke triangles and the Wigner ansatz for this quantity (thick solid hues). Top - 
the first 10000 levels for the triangle with n = 5, middle - the first 20000 for the 
triangle with n = 7, bottom - the first 10000 levels for the triangle with n = 12. 
Dashded line - the same quantities but for the triangles with angles lOvr/llQ, 
lOvr/Tl, 207r/99, respectively. Thick solid lines at each graph are the difference 
between the true GOE prediction and the Wigner ansatz. 



the difference between the true GOE formula and the Wigner ansatz. From 
the figure it is clearly seen that spectral statistics for the Hecke triangles 
is quite close to the conjectured statistics of the Gaussian Orthogonal En- 
semble (GOE) of random matrices. To have an estimate of statistical errors 
we compute numerically the same quantities (dashed lines at Fig. IT3|) for 
non-tesselating triangles of the Hecke type which have two angles 7t/2 and 
but instead of the angle tt/tt, as for the true Hecke triangle we take a certain 
angle 7„ sufficiently close to it. For n = 5, 7, 12 we choose respectively 
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Notice that the difference between 7„ and vr/n is quite small 
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For all cases (except possible small deviations for n = 12 which has the largest 
multiplicities) the nearest-neighbor distributions for the Hecke triangles agree 
well with the curves for non-tesselating triangles 

These (and others) figures demonstrate that the spectral statistics of the 
non- arithmetic Hecke triangles (even with quite large degeneracies of periodic 
orbit lengths) at small distances is rather well described by standard random 
matrix ensembles. 

The contradiction between the observed random matrix statistics of the 
Hecke triangles and deviations of correlation functions due to large multiplic- 
ities of periodic orbits (cf. (IHUj) was partially resolved in jlj. In this paper 
it was demonstrated that the diagonal approximation can, strictly speaking, 
be applied only for very small values of t < ti. If the mean multiplicity 
increases like g{l) ~ e'^' with a certain constant A < 1/2 from |3] it follows 
that the time of applicability of the diagonal approximation has the following 
estimate 

1 In A; , , 

During this time the form factor increases exponentially but it can reach only 
a value of the order of 

K{ti) ~ A;-(1-2A)/{1-A) _ (84) 

For arithmetic systems A = 1/2 (see 0) and the form factor for the time of 
applicability of the diagonal approximation becomes of the order of 1 which 
explains the Poisson character of their spectral statistics. But for all non- 
arithmetic groups A is less than 1/2 and the form factor in the diagonal 
approximation increases only by a negative power of k. Therefore in the 
semiclassical limit k oo there is no apparent contradiction between ob- 
served GOE-type local statistics and the change of correlation functions due 
to large multiplicities of periodic orbits. 

These arguments suggest that two-point form factors for non-arithmetic 
Hecke triangles has the form indicated at Fig. EI The peak at small values 
of t is due to large multiplicities of periodic orbits. The magnitude of this 
peak and its position seem to decrease for large k. 

Though deviations from standard statistics should be small when k — >■ 
oo the peak indicated at Fig. HH may influence the large distance spectral 
properties like the number variance (see e.g. f3^). At Fig. I15IIT7I we present 
the number variance for the Hecke triangles with n = 5, 7, 12 together with 
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Figure 14: Schematic form of the two-point correlation form factor for non- 
arithmetic Hecke triangles. Dashed line is the continuation of the GOE form 
factor to small values of t. 

the corresponding values for non-tesselating triangles (jHT|) . Due to large 
statistical errors in the computation of the number variance S^(L) we found 
convenient to plot at the figures not itself but its averaged value defined 

in the following way 



To demonstrate the evolution of the number variance with increasing the 
energy at all figures we present pictures for the averaged number variance 
with different number of levels. 

For all non-tesselating (generic) triangles the number variance follows the 
the GOE prediction for small values of L and then saturates, as it should be 
for dynamical systems PP, j2]. For the Hecke triangles the number variance 
at small L also agrees with the GOE formula but then it becomes bigger than 
this reference expression and only later it saturates but at a value different 
from the one of the corresponding very close-by non-tesselating triangle. 

This overshooting looks as a direct confirmation of the conjectured form 
of the two-point correlation form factor (cf. Fig.ll4p but careful calculation of 




(85) 
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Figure 15: Smoothed number variance for the Hecke triangle with n = 5 (top 
graphs). Different type of hues corresponds to the number variance computed 
from the first 2000 • k levels. Solid line - A: = 1, dotted line - k = 2, dashed line - 
k = 3, long dashed line - k = 4, dot-dashed line - A; = 5. Thick line is the GOE 
prediction. Bottom graphs represent the same quantities for the non-tcsselating 
triangle with angle 75 = 207r/99. For clarity top graphs are shifted up by .2 unit. 
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Figure 16: The same as at Fig. El but for the Hecke triable with n = 7 (top 
graphs). Different types of lines correspond to the averaged number variance 
computed from the first 4000 • k levels. Bottom graphs are calculated for the 
non-tesselating triangle with angle 77 = lOvr/Tl. 
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Figure 17: The same as at Fig. EJbut for the Hecke triangle with n = 12 (top 
graphs). Bottom graphs are calculated for the non-tesselating triangle with angle 
712 = lOvr/lig. 
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this quantity requires a ressumation of, at least certain, non-diagonal terms 
and is beyond the scope of this paper. 

7 Summary 

We demonstrate both numerically and analytically that, at least, certain 
non-arithmetic Hecke triangle groups have exponentially large multiplicities 
of periodic orbit lengths. 

In groups under consideration matrix elements of group matrices are in- 
tegers of an algebraic field of a finite degree and each group matrix gives rise 
naturally to q different lengths corresponding to q different isomorphisms of 
the basis field. The main ingredient of our approach to the problem of peri- 
odic orbit length multiplicity is the investigation of the joint distribution of 
periodic orbits with all q transformed lengths fixed. 

We conjecture that this distribution has a scaling form PU)) and find the 
scaling exponent numerically. For Hecke groups (j24p with only one non- 
trivial isomorphism the general inequality is sufficient to demonstrate 
an exponential increase of the multiplicities. Multiplicities obtained by this 
method are in a good agreement with direct numerical calculations. 

For general Hecke triangle groups we are not aware of analytical condi- 
tions of the existence of necessary solutions. In all investigated cases (except 
for the groups (j24p with only one non-trivial isomorphism) the increase of 
multiplicities numerically is too small to be observed from direct calculations 
of periodic orbits but the data fluctuate around a value bigger than 2. 

The spectral statistics of non-arithmetic Hecke triangles agrees with the 
GOE statistics at small distances but deviates from usual expectations at 
large distances. 
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Appendix A 

To prove the inequality ()66p it is slightly more convenient to describe conju- 
gacy classes of the Hecke triangle group matrices not by the code discussed 
in Section |21 but by a code proposed in |Tn] better suitable for analytical 
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calculations. In this code the letters for the orientation preserving subgroup 
of Gn are the following matrices 



gk{an) = U''-'T , k^l,...,n-l, (86) 

where U = TS and matrices T, S are the translation and inversion matrices 
which generate the whole group Gn 





with a = 2 cosvr/n. 

As in the previous code periodic orbits for the Hecke group G„ (with unit 
determinant) are free words of letters gk, the only restriction being that all 
cyclic permutations of a word give one orbit. 

It is easy to check (e.g. by induction) that 

9k=(''' ""rO (88) 



where = ajk(2cos^) are the values of the Chebyshev polynomials of the 
second kind 

sin(A;^) , , 

afe(2cos^) = — ^ . (89) 

computed at 6 = ir/n. 

Let us introduce the following definition. We say that a function f{x) 
has the H-property with a separating point h if 

\f{x)\ < f{h) for all \x\<h. (90) 

The importance of this notion follows from the fact that if /i [x) and /2 {x) 
both have the H-property with separating point h then fi{x)f2{x) and fi{x) + 
/2(x) also have the H-property with the same separating point. In particular, 
if one has a set of matrices whose elements all have the H-property with 
a separating point h then all products of these matrices also have the H- 
property with the same separating point. 

Let us prove first that matrix elements Ofe with k — 0, . . . ,n have the 
H-property with the separating point h — 2cos(7r/2n). Indeed 

|afe(2cos^)| < ^ . (91) 
sm9 
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When < ^ < 7r/2 the equahty sign in this inequahty holds at the points 

Om = (92) 

where m is an odd integer 1 < m < k. 
Therefore when 6 > tt /2k 

|afc(2cose)| < ^ <afc(2cos;^) . (93) 

Sill U Z,rC 

Moreover, aA;(2cos6') is a decreasing function when < 6 < 7r/2A; and 
afc(2cos6') < afc(2 cos7r/2n) when n/2n < 9 < Tr/2k (and, of course, k <n). 
Together these two inequahties prove that with k <n have the H-property 
with separating point h = 2cos(7r/2n). 

Second, 0^(2 cos vrp/ra) = a„_fc(2cos7rp/n) for odd integer p. It means 
that for all isomorphisms (j^^ ak{x) equals an-k{x) and only with k < n/2 
are independent. Hence, for all isomorphisms of the defining equation 
can be considered as polynomials of degree not greater than n/2 and one 
can choose for all 0^(2 cos vrp/n) with k < n the same separating constant 
h = 2 cos(7r/n). 

Third, as was stated above, all matrix elements obtained by taking the 
products of arbitrary number of matrices gt also have the H-property with 
the separating constant h = 2cos{'K/n). 

Combining all these arguments one proves that for all isomorphisms of 
the basis field traces of the Hecke group matrices have the H-property with 
the same separating constant. Because 



nk 
cos — 

n 



71" 

< COS - (94) 

n 



for all k 7^ 0,1, tt, one gets that modulus of traces of the Hecke triangle 
group matrices decrease for all non-trivial isomorphisms of the basis field 
thus proving the inequality ()66|) . For matrices with determinant equal — 1 the 
same inequality follows by computing the square of such matrices because 
the periodic orbit length for the square of any matrix is twice the length 
corresponding to the initial matrix. 

After this paper has been completed we become aware of Ref. [TS| where 
the inequality (jUUI) was proved for all groups which permit the so-called mod- 
ular embedding. From Ref. [7j it follows that all triangle discrete groups 
belong to this class. Therefore, the inequality ()66|) is valid for all triangle 
groups (and not only for the Hecke triangle groups considered in this paper). 
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Appendix B 



The purpose of this Appendix is to give arguments in favor of the represen- 
tation (jHHj) of the joint probabihty density of periodic orbits with all trans- 
formed lengths fixed. 

For discrete groups periodic orbits can be obtained from product of cer- 
tain number of matrices. Let us consider in a given code the product of n 
basis matrices 

A{n) = An-An-i---A^. (95) 
The total number of matrices with n symbols for a general code is exponential 

phn 

p(^) ^ (96) 
n 

where /i is a constant called the topological entropy. 

The length of periodic orbit is related with matrix A asymptotically as 

/ = 2 In Tr A . (97) 

Therefore, matrices representing periodic orbits can be considered as the 
result of a random process where matrices Ak are chosen randomly from a 
code grammar. The probability distribution of lengths for products of n 
such matrices is defined as the ratio of number of matrices with lengths in 
the interval / + dl] divided by total number of matrices. 

This distribution under quite general conditions jH] , [ID] has the Gaussian 
form 

P^{1) = -^^e~^^ (98) 



where 

In = Aon , al = aln . (99) 

One of possible applications of such distribution is the calculation of the total 
density of periodic orbits of length / (see e.g. |T11) 

/phn 
—Pn{l)dn . (100) 
n 

When / — ^ oo the integral can be computed in the saddle point approximation 
and the total density of periodic orbits has exponential asymptotics 

p{l) = - (101) 



33 



where 

K = ^ ° . (102) 

For groups considered in the paper all matrix elements belong to an algebraic 
field of a finite degree which has q different isomorphisms. It means, in 
particular, that each product of n group matrices A{n) as in (j95p gives rise 
to q different lengths li{n) obtained by applying each isomorphism ipi to A{n) 

k{n) = 2\n\i:i^i{A{n))\ . (103) 

On the other hand ipi{A{n)) can be obtained as the product of n transformed 
matrices ipi{Ak) as in Therefore according to the above theorem each 

variable U is a random variable whose distribution also has asymptotically 
the Gaussian form 

= -JK^^ (104) 
V Zirnai 

with certain parameters Aj and cxj having the meaning of the mean value and 
the variance of /j(n). 

Let us conjectured that the mutual distribution of all li{n) together is 
also Gaussian 



^"^^"^^ = ("i t.M^,{h - X^n){l, - A,n)j (105) 

with a certain positive definite matrix Mij. 

Analogously to Eq. (jlOOj) the total density of orbits with fixed k is 

roc phn _^ 

p{l)= —Pnil )dn . (106) 
Ji n 

As above this integral can be computed in the saddle point approximation 
valid at large / and the result is 



= (2.)(.-w"- 2k))U^ "MB- ^MC - 2h) ) (107) 



where 



A=J2 Mijhlj , B = J2 M^.kXj , C=J2 MijXiX, . (108) 

i,j=l i,j=l i,j=l 
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The exponent in (|1(J7|) is an homogeneous function of Ij and after the rescahng 
Xj = Ij/h one obtains the scahng ansatz (j46p with a specific scahng function 
which leads to exponential asymptotics of the joint probability distribution 
as it seems suggested by numerics (cf. Figs. |S]|H1). 

The main drawback of such approach is that the theorem about the Gaus- 
sian form of the distribution of the product of n random matrices is valid 
only near the maximum of the distribution. But the term e^"' in (ll(J(J|) and 
fll07p shifts the saddle point far from the maximum and there exist no general 
arguments which would imply the smallness of corrections to the parabolic 
form of the exponent. For certain groups and special codes it seems that one 
can ignore such corrections in a region of interest but, in general, corrections 
are large and one has to rely on the numerics as it was done in the main part 
of the paper. 
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